Collective oscillations of a Fermi gas near a Feshbach resonance 
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A sum rule approach is used to calculate the zero temperature oscillation frequencies of a two 
component trapped atomic Fermi gas in the BCS-Bose Einstein condensation crossover region. These 
sum rules are evaluated using a local density approximation which explicitly includes Feshbach 
molecules. Breathing modes show non-monotonic behavior as a function of the interaction strength, 
while quadrupole modes are insensitive to interactions for both spherically symmetric and axially 
symmetric traps. Quantitative agreement is found with experiments on atomic 6 Li system and with 
other theoretical approaches. 



I. INTRODUCTION 

The successful trapping of fermionic atomic gases in 
magnetic and optical traps and their subsequent cooling 
to quantum degeneracy has attracted much interest. Ex- 
perimentalists wield large degrees of experimental control 
over these low temperature and low density gases, includ- 
ing the ability to manipulate the effective atom-atom in- 
teraction by tuning an external magnetic field 0,00,0- 
The Zeeman shift from the magnetic field adjusts the en- 
ergy of a two atom bound state relative to the scatter- 
ing continuum. Around the Feshbach (FB) resonance, 
where the binding energy is exactly equal to the energy 
of colliding atoms, a small change in magnetic field can 
have a dramatic effect on the effective interaction 0. 
On one side of the resonance (above the resonance) the 
atomic gas has attractive interactions, while on the other 
side (below the resonance) effective interactions are re- 
pulsive. At resonance, the scattering length of the atoms 
(a) is divergent and one expects universal behavior 0. If 
the temperature is low enough, then fermionic atoms can 
pair and undergo Bardeen, Cooper and Schrieffer (BCS) 
condensation on the attractive side (BCS regime) of the 
resonance. On the repulsive side (BEC regime), there ex- 
ists Bose molecules (FB molecules) which are short-range 
fermionic atom pairs and these molecules can undergo 
Bose-Einstein condensation (BEC) at low temperature. 
By sweeping the external magnetic field through the FB 
resonance, one expects a continuous crossover between 
these limits @ H G3 GH H El 

Here we study the collective oscillations of a Fermi gas 
in the crossover region. As a readily accessible dynam- 
ical observable, collective oscillations of a trapped gas 
can yield crucial information on the state of the system. 
Recently, both transverse and axial breathing modes of 
a superfluid 6 Li atomic system have been measured in 
a highly anisotropic trap 0, 0, 0, 0. Several the- 



oretical attempts have been made to understand these 
collective oscillations [IlllSiJIIllIllIlilllil. 

Stringari used a hydrodynamic theory to calculate 
the collective oscillations far from resonance (|a| <C 1) 
and at unitarity (|a| — > oo). By interpolating from these 
results he predicted that the frequency of the transverse 
breathing mode in a highly elongated trap should exhibit 
a non-trivial dependence on the scattering length jl9j . 
This prediction was further supported by studies based 
on semi-empirical forms of the equation of state, time 
dependent density functional approaches, and asymp- 
totic expansions 0, |U |H Hi [H HI . Additionally, 
with explicit inclusion of Feshbach molecules, Ohashi 
et al |2y| used a Hartree-Fock-Bogoliubov (HFB) ap- 
proach to study the collective oscillation in the BCS-BEC 
crossover region. While their theory predicts correct re- 
sults in the weak coupling BCS limit, their results in the 
weakly repulsive BEC limit are those of a non-interacting 
(rather than interacting) Bose condensate. 

Unlike some of these previous studies, our approach, 
based upon the "compressibility sum rule" and the "f- 
sum rule", describes this system in all regions of the 
phase diagram. It is more flexible than the methods 
based upon isotropic indices, explicitly includes Feshbach 
molecules, and correctly predicts all limiting cases. 

The paper is organized as follows. In section II, the 
two channel coupled fermion-boson model is introduced. 
Following the pioneering work in the context of supercon- 
ductivity by Eagles , Leggett , Nozieres and Schmitt- 
Rink we simplify this model by using a mean field 
theory described in section III. In section IV, we use 
sum rules to derive the collective breathing modes and 
quadrupole modes for both spherically symmetric and 
axially symmetric traps. Finally, we discuss our results 
with comparison to the experimental and other theoret- 
ical predictions in section V. 



II. THE MODEL 
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We consider a dilute gas of interacting Fermi atoms 
in two hyperfine states trapped in an external poten- 
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tial V(f). Since the atomic gas is cold and very di- 
lute, the inter-atomic distance is large compared to the 
range of the atomic potential and s-wave scattering dom- 
inates over all other scattering channels. Therefore, we 
neglect interactions between fermions in the same hyper- 
fine states and approximate the attractive inter-atomic 
potential between fermions in different hyperfine states 
by a contact potential — Uba,i C S(rf — r). Furthermore, we 
take into account coupling between molecular bosons as- 
sociated with the Feshbach resonance and Fermi atoms 
by a local conversion term, proportional to <?barc- The 
Feshbach resonance is achieved by tuning the thresh- 
old energy 2^b a rc which varies linearly with magnetic 
field. The atomic gas is then described by the Hamilto- 
nian El El El E5 



ft 2 V 2 
2m 



■fM) + V(f) Vv(f) 



+ J d 3 r<f>l(r) ( - + 2^ barc - 2 Mo + 2V(r)^j 0(f) 

+9b„» J d 3 f^(f)Mm'1 If) + h.cj .11) 

where m(M — 2m) is the mass of the Fermi atoms (Bosc 
molecules), and (M)(fJ>M = 2/xo) is the Fermi (Bose) chem- 
ical potential. The field operators -(/v(f)[0(f)] obey the 
usual fermionic (bosonic) anticommutation (commuta- 
tion) rules, and describe the annihilation of a fermion 
(boson) at position fin the hyperfine state a = \ or J,. In 
general, the trapping potential for Fermi atoms is taken 

as V(r) =Yl a 3 — T^ r aTf3. For Bose molecules, the trap- 
ping potential is 2V(r). Notice, the trapping frequen- 
cies uJa/3, a, P = x,y,z are the same for both species 
of fermions and molecular bosons. In this paper, we 
consider only spherical symmetric and axially symmet- 
ric traps. For the case of a spherically symmetric trap, 
uj a p = Lo5 a jj 1 where S is the Kronecker delta. For axi- 
ally symmetric traps, aj a /3 = Lu±5 at p for a, f3 — x,y and 
= bJ z 5 at p for a,f3 = z. 

III. MEAN FIELD EQUATIONS 

The gap equation and the number equation are ob- 
tained by using a mean field approximation. Introducing 
two local mean fields for the BCS condensate and BEC 
condensate respectively, A(f) = u^e^i^-ip^r)) and 
4>m(r) = (0(f)), the mean field Hamiltonian is written as 



Hmf — 



ft 2 V 2 
2m 



fj,(f) )4> a {r) 



d 3 f0 f (f) 



2M 



2 M (f) U(f) 



h.c 



(2) 



where we defined a local composite order parameter 

A(f) = A(f) - 5b are<Mf) = <L W T (0) with 

<tc(r) = Ubaxe ~ 2 M (^-2l barc ■ Where ( A ) represents 
the expectation value of the operator A. The local chem- 
ical potential fi(r) = fi — V(r) is to be treated with local 
density approximation (LDA). Within the LDA, we as- 
sume that the density varies so slowly so that it can be 
treated as locally uniform. Therefore, we determine the 
solutions of the equations at each point in space as for a 
homogeneous system. This approximation is reasonable 
for the current experiments as the number of atoms in 
the trap is very large [2sl Efll ] . Also, we assume that at 
zero temperature all the molecules are Bose-condensed. 
The two order parameters introduced before are not in- 
dependent and they are strongly coupled to each other 
by 



ffbarcA 
^barc 



+ (2tWe - 2fji)<t> m = 0. 



(3) 



We have suppressed the position dependence of the local 
quantities in this equation and the following. Notice, as 
a result of this strong coupling between A(f) and m (f), 
both BCS and BEC condensates are non-zero through 
out the entire BCS-BEC crossover region. 

Diagonalizing the mean field Hamiltonian by a Bogoli- 
ubov transformation, the usual BCS gap equation and 
number equation at zero temperature are given by 



1 = u 



bare 



n F (r) 



(r) 



d 3 k 
J2nf 



d 3 k 1 
(2tt)3 2E k 



(1 



E k 1 



Here € k = h 2 k 2 /2m is the kinetic energy, E k 



(4) 



(5) 



(efc — m) 2 + A 2 is the BCS-Bogoliubov excitation en- 
ergy, and UF^r) is the density of the Fermi atoms. The 
total density is the sum of fermionic atom density and 
Bose molecular density, n(r) = np(r) + 20 m (f) 2 . Be- 
cause of the short range nature of the interaction, the 
gap equation shows an ultra-violet di vergenc e. Using a 
standard regularization procedure |l(lllll,ll4Ll3C| to sub- 
tract off the divergences, a renormalized effective inter- 
action Upff(f) is introduced through the relation — = 

otf 1 , = — f ,i % o~ ■ Then the renormalized BCS gap 
equation is written as 
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1 = Ueff(f) 



d 3 k 1 J_, 
(2tt) 3 ( 2£ fe ~ 2e k ' 



(6) 



where renormalized effective interaction is written in 
terms of renormalized pairing interaction u, renormal- 
ized FB resonant strength g and renormalized detuning 



v as u cS (r) 



Notice, we 



1 2[i{r)~2v m 

have introduced a position dependent scattering length 
a(r) which crucially depends on the local chemical po- 
tential /x(r) and the threshold energy 2v (See FIG. 1 be- 
low). The composite order parameter A(r) is now writ- 
ten in terms of renormalized g as A(r) = A(r) — g<j) m {r). 
The renormalized non-resonant pairing interaction u and 
renormalized FB resonant strength g must be extracted 
from experiments. Assuming the atom scattering length 
near FB resonance takes the form a s = ab(l — B W Bg ), 
with aj, the background scattering length, B , the reso- 
nant magnetic field, W, the resonance width, they are 

given by u = ^ and g = y ^MM wnere w is 
the difference of the atomic magnetic moments between 
the closed and the open scattering channels. From exper- 



exper- 
el E3, 



iments on Li atomic system by Bartenstein at 
we see that it ~ 9.0 x 10 7 fcs nm 3 /iK and g ~ 1.7 x 10 6 fc_B 
nm 3,/2 //K, corresponding to the broad resonance of e Li 
system at magnetic field Bo ~ 837G. In addition to this 
broad resonance, this system shows a narrow resonance 
at B ~ 545G which gives g ~ 2.4 x 10 3 fc B nm 3 / 2 ^. 
For the broad resonance, the crossover occurs at \u\ ^> e/, 
while at the narrow resonance, it occurs at |z/| ~ e/, 
n 2 k 2 

where £/ = 2m / is the Fermi energy. However, when 
viewed as a function of scattering length, one finds the 
same crossover physics irrespective of the width of the 
resonance 33]. In this paper, we have done our calcula- 
tion for the broad resonance and present our results as a 
function of the dimensionless parameter (fc/a s ) _1 , where 
a s is the scattering length at the center of the trap de- 

2 Ah 2 

fined through the relation, u — „ 9 „ = — a ° , hence 
our results must be applicable to the narrow resonance 
as well. 

Eq. (5) and Eq. (6) are to be solved with the constraint 
that the total number of atomic particles in the trap is 
N = J <firn{r). The size of the cloud tq is determined 
by the condition n(f — fb) = 0. Depending on the sign 
of the scattering length, this condition gives two different 
formulas for determining the size of the cloud. For a(r) < 
(BCS regime), the boundary is given by /i(r = r*o) = 
and, for a(f) > (BEC regime), it is given by /z(f = 
r*o) = — m a[ r - a yi ■ I n the weak coupling limit (fcja — > 0~ 
or equivalently v — > +oo), Eq. (5) and Eq. (6) reduce 
to the standard BCS results, n(r) cx /i 3 / 2 and y, <~ e/. 
When kfa — * + (or equivalently v — oo), they lead to 



a Bose gas with n(r) cx; |/i| 2 (- 



T) and /x 



'a^/2m|/i| ; 2ma 2 ■ 

At unitarity (a = ±oo), physical quantities of the system 



must be independent of the scattering length. Therefore, 
the system shows universal behavior with n(r) cx /z 3 / 2 [|J. 

We plot the spatial variation of the inverse scattering 
length and density in FIG. 1 in a spherical trap, using 
parameters from the two resonance in e Li. We choose 
v ~ 4e/ for the narrow resonance and v ~ 10 4 e/ for the 
broad resonance so that a(0)fc/ is the same in each case. 
The large variation in et(r) for the narrow resonance is 
due to the fact that v ~ e/. 



re 




FIG. 1: Spatially dependent inverse scattering length 
{fc/a(r)} _1 for spherically symmetric trap as a function of 
r/lo, where Iq and kj are oscillatory length and Fermi wave 
vector at the center of the trap respectively. Dotted line rep- 
resents the broad resonance while solid line represents the 
narrow resonance. Inset shows the scaled density variation at 
the same parameter regime. 



IV. SUM RULES AND THE COLLECTIVE 
OSCILLATIONS 

We use a sum rule approach to study the propagation 
of collective modes in harmonically trapped Fermi gas 
in the presence of FB resonance. Sum rules are derived 
from the linear response theory and allow us to compute 
rigorous upper bounds to the energies of the collective 
oscillations at zero temperature. The response function 
x(z) at complex frequency z is written in terms of the 
imaginary part of the response function x' 1 as 



X(z) = 



7T (iO — Z) 



(7) 



The explicit determination of the dissipative component 
of the response function x"( w ) requires the full solution of 
the Hamiltonian. However, useful information on the be- 
havior of the system can be derived by using the method 
of sum rules. Expanding the response function in terms 
of 1/z, moment expansion of x( z ) is given by 



4 



where the pth moment is defined as 



(8) 



lp = J ^V x ». 



The collective oscillation of the system is then bounded 
by \J 2 l^l- With the sum rule approach, it is 
possible to compute these moments without evaluating 
response function. Considering quadrupole excitations 
Qap = f d 3 rr a rp(^ l(T 4>t' l l J cr + 2<^</>) and the perturbed 
Hamiltonian H pert = JJ2 a /3 d 3 r\ a pr a rp(^ a iptipa + 
2</J<f>), the response function is written as x"s ~ 



where A, 



a/3 



2 mUJ al 3 



dX-, s ' 

and a, (3 run over the Cartesian 

indeces x, y and z. The moments, (rn_i)"^ and 
are given by 



(m-i)# = 
= 0) = 



1 d(Q a p) 

muj 1 s duj-ys 



(10) 



and 



Mj = / ^Im(^) = {[[Q af3 ,H},Q^}) 
— (5 ai (Qp5) + 6 a s(Qp 7 ) + S/3 S (Q ai ) + 5^{Q a s))(ll) 



where [A, B] represents the commutator between the op- 
erators A and B, The squares of the collective mode 
frequencies are approximated by the eigenvalues of the 
matrix M = (rn,—i) (mi) ■ 



m 



4 (r 2 



3 2<£l 

1 



d(x 2 
db v 



±fd 3 r(x 2 



2 \2 9n(r) 



(13) 



where, (r 2 ) = J d 3 rr 2 n(f) is the expectation value of 

(9) 

the operator r 2 with the total atomic density n(r). In 
the weak coupling BCS limit, the number of molecu- 
lar bosons can be neglected, hence one can show that 
Vb and (r 2 ) cx As a result, the breathing 

mode and the quadrupole mode approach to the expected 
results to m — > 2uj and uj q — > ^/2us as in the case of a 
non interacting Fermi gas. In the opposite limit, the 
number of Fermi atoms can be neglected, so that one 

can show that u„ ( — ,I l — 1)§ oc hi and 

'r.2\ ~- "- 1 6~ 1 ( 2to ^ 2 — 1) which resulted in w m — > \/jkd 



CX (Ltot 



and uj q — > \p2~Lo. Notice, in this limit, the mode frequen- 
cies are identical to the result obtained by the Thomas- 
Fermi limit of a Bose gas j3^- As opposed to the HFB 
theory |2fij . our approach retains effective molecular in- 
teraction at the weak repulsive BEC limit and recovers 
hydrodynamic results for both limits j^HH^- Also, at 
unitarity, the chemical potential dependence on the den- 
sity is the same as that of a non interacting Fermi gas, 
hence the mode frequencies are expected to be the same 
as those of non interacting limit of the Fermi gas. We 
numerically evaluate Eq. (12) and Eq. (13) in the LDA. 
The breathing mode frequencies in the crossover region 
are plotted in FIG. 2 as a function of dimcnsionless pa- 
rameter (fc/a s ) _1 . In addition to finding the expected 
results as (fcya s ) _1 — ► ±oo, we find a dip near unitarity 
on the BCS side. 

Our numerical calculation confirms that the 
quadrupole mode calculated from Eq. (13) is inde- 
pendent of the interaction and the aspect ratio of the 
trap (see below), and its frequency is y/2u> through out 
the entire BCS-BEC crossover. 



Spherically symmetric trap 



Axially symmetric trap 



In the case of a spherically symmetric trap, V(r) = 
^mu> 2 r 2 — br 2 . Solving for the eigenvalues of the matrix 
M, the breathing mode uj m , and the quadrupole mode u> q 
associated with eigenvectors {1, 1, 1} and { — 1, 1, 0} are, 
respectively, given by 



4 (r 2 ) 
m a < r2 > 

db 



(12) 



and 



In the case of axially symmetric trap, V(f) = 
imu^(p 2 + ^ z2 ) = bj_p 2 + b z z 2 , where p 2 = x 2 + y 2 . 
For this case, solving for the eigenvalues of the matrix 
M, the breathing modes u± are given by 



1 



m(Q 1 Q 2 -2Q 3 )\2 

ha 2 ^ 

•4 

1 



t;(p 2 )Q2 + (z 2 )Qi 

2 . 



±{- a (p 2 ) Qi + (z 2 ) Ql 



+ ^(p 2 )(z 2 )(SQ 3 -2Q 1 Q 2 )}- 



(14) 
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FIG. 2: Zero temperature breathing mode frequency of a two 
component Fermi gas in a spherically symmetric trap within 
the BCS-BEC crossover region. The frequency is given as a 
function of a dimensionless parameter (fc/a s ) _1 . Our the- 
ory predicts uj m /u — > 2 at the weak coupling BCS limit 
[(fc/Os) -1 — > — oo] and w m /cu — ► v5 at the weakly repul- 
sive BEC limit [(fc/a a ) _1 — > +oo] as one expects from hy- 
drodynamic theory. At unitarity, uim/u -> 2 as required by 
universality. 



where we define, Qi — 2~WT : ® 2 = ~~db7~ an< ^ = 
I " ^dbj • Q ua drupole mode w ? is given by 



the case of highly elongated traps, u>+ — > 2a; j_ and 
cj_ — ► y^|w z . At unitarity, collective frequencies are ex- 
pected to be the same as those at the weak coupling BCS 
limit. The mode frequencies uj + and w_ in the crossover 
region are presented in FIG. 3 and FIG. 4 for the case of 
A = 0.001. The behavior of the modes are the same as 
that of a spherically symmetric trap, agreeing with the 
hydrodynamic theory of a non interacting Fermi gas and 
a large N Bose gas at the asymptotic limits, while dis- 
playing a dip near the unitary limit on BCS side. A sim- 
ilar behavior of the collective oscillations has been recov- 
ered using either a scaling ansatz based on a mean field 
description of the BCS-BEC crossover or the gen- 
eralized Hylleraas-Undheim method ji^. As suggested 
by Hu et al [2(|, these dips in the collective frequen- 
cies may be attributed to the fact that pairing enhances 
the compressibility of the atomic gas. The measured ax- 
ial breathing modes have quantitative agreements in the 
BEC side with experiments carried out by Bartenstein et 
al |l6j . The measured transverse breathing modes are in 
quantitative agreements with the experiment carried out 
by Kinast et al [l5| . For comparison, we plotted exper- 
imental results from Ref. ^3 and Ref. 0] as solid 
circles in FIG. 3 and FIG. 4. In section V, these experi- 
mental results are compared with our prediction in more 
details. 

As in the case of a spherically symmetric trap, the 
quadrupole mode ui q — > y/2ui± and our calculation con- 
firms that it is independent of the interaction. 



4 (p 2 



y 



2\2 Sn(r) 



(15) 



The expectation values (....) are defined as before. In 

the weak coupling BCS limit, the trapping frequencies 

i 

dependence on the chemical potential as po oc b±_b z . At 
this limit, the frequencies dependence on the expectation 

values of the operators p 2 and z 2 are (p 2 ) oc b z 2 

Then Eq. (13) yields w| 

4c^(^ + g±^ + £_3>). For highly the elon- 

gated traps (A <C 1) in which most of the experiments 



and (z 2 ) cx pob ± b z 2 ■ 

,2/5 i A _i_ /~25 , A 2 2A 



luj_ and 10- 



In the 



are carried out, u> + - 

opposite limit (weakly repulsive BEC limit), the trapping 
frequencies and the chemical potential are related by the 



expression, \po\ 2 (- 



1)( 



h 

2ma' 2 |/x 1 



1)9 oc bxb'i 



In this weakly repulsive BEC limit, the expectation val- 
ues of the operators p 2 and z 2 , the chemical poten- 
tial and, the trapping frequencies are related as (p 2 ) cx 



ay/2m\fj, \ 



ImoK 

then yields u>± 



1)( 



2ma 1/j.q I 



b^bl 



-1)( 



2ma' 2 I 



E, ,2 / 2 1 6a _i_ I 4 1 yA^ 4A 
^lU + 10 ± V 25 + TOO 25 



and (z 2 ) oc 
TheEq. (14) 
For 



V. DISCUSSION AND CONCLUSIONS 

Our results are in quantitative agreement with the ex- 
periments carried out by Kinast et al (ljj and Barten- 
stein et al [l(j on the Fermi gas of 6 Li atoms. While 
transverse breathing modes have been measured in both 
experiments, axial breathing modes have been measured 
only by Bartenstein et al |16| . The measured finite tem- 
perature axial breathing modes by Bartenstein et al 0] 
initially show a decrease in the collective excitation fre- 
quency as the magnetic field is swept from the weakly 
repulsive BEC limit towards the BCS limit. Experimen- 
tally a dip is found at (fc/as)^ 1 ~ —0.5 after which 
there is a striking increase in the oscillation frequency 
(see FIG. 3). The location of the dip is in excellent 
agreement with our prediction which shows a similar 
dip at (kfas)^ 1 ~ —0.46. The measurements at the 
weakly repulsive BEC limit (u>- = \Jb/2uu z ) and at 
unitarity (cj_ = y/l2/5u> z ) are in quantitative agree- 
ment with our prediction. However, the measured collec- 
tive axial breathing mode frequencies at the weak cou- 
pling BCS limit approach 2uj z , while, as with most other 
theories, we predict that oj_ — > y/T2/5u> z . In fact at 
(kfas)^ 1 r~j —0.67, the upper bound calculated from our 
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1.565 



1.555 



1.545 



1.535 




-1.5 , 1.5 

FIG. 3: Zero temperature axial breathing mode frequency of 
a two component Fermi gas in a highly elongated (A = 0.001) 
trap within the BCS-BEC crossover region. The frequency 
is given as a function of dimensionless parameter (fc/a s ) _1 . 
Notice, at the weak coupling BCS limit [(fc/a s ) _1 — > -co], 
lu~/lo z — > y^f and at the weakly repulsive BEC limit 
[(fe/a s ) _1 — > oo], u-/lo z — » v/¥ as one expects from hydro- 
dynamic theory. At unitarity, tu-/tu z — > y^j? as required by 
universality. The solid circles are finite temperature experi- 
mental results from Ref |T3| . 

sum rules begins to be below the measured axial mode 
frequencies. As suggested by Bartenstein et al ^(J, this 
jump may indicate a breakdown of the hydrodynamic 
theory, and is undoubtedly related to pair breaking as 
at this point the mode frequency is comparable to twice 
the gap. If the experiment is done in the ranges of densi- 
ties and temperatures where the threshold energy of the 
two-particle continuum 2E g (where E g = A for /i > 

and E g = ■J fi 2 + A 2 for \i < 0, is the single particle 
excitation gap) is smaller than the collective energies, 
then the collective modes can become overdamped, leav- 
ing the system to simply oscillate at a multiple of the 
trap frequency. Since we work in LDA and make mean 
field approximations our upper bound is not rigorous in 
this region. 

As shown in FIG. 4, our prediction for the transverse 
breathing mode frequencies are in quantitative agreement 
with the finite temperature experimental measurements 
by Kinast et al |l5| . The experimental data points at 
(fc/a s ) _1 ~ —0.42 represent different temperatures. Sur- 
prisingly, we have a good agreement with higher tem- 
peratures at (fc/a s ) _1 ~ —0.42. Our predictions do not, 
however, agree with the transverse breathing mode data 
reported in Ref. [lfij . which show an abrupt change and 
large damping close to the unitary limit. This disagree- 
ment may be due to ellipticity in their optical trap |37j| . 

Finally, our results can be compared to the other the- 
oretical approaches. We have very good agreements in 
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FIG. 4: Zero temperature transverse breathing mode fre- 
quency of a two component Fermi gas in a highly elon- 
gated (A = 0.001) trap within the BCS-BEC crossover re- 
gion. The frequency is given as a function of dimensionless 
parameter (fc/a s ) _1 . Notice, at the weak coupling BCS limit 
[(fc/ds) -1 — > — oo], cj+/lu± — > \/t and at the weakly repul- 
sive BEC limit [(fc/a s ) _1 — > oo], lo+/uj± — > 2 as one expects 
from hydrodynamic theory. At unitarity, lj+/uj± — > as 
required by universality. The solid circles are finite temper- 
ature experimental results from Ref [l5f . The multitude of 
experimental data points at (fc/a^)^ 1 ~ —0.42 represent dif- 
ferent temperatures. 



the entire BCS-BEC crossover region with the theories 
proposed by Hu et al |2(J using a scaling ansatz together 
with a mean field description of the BCS-BEC crossover 
and Kim et al |22j using an approach based on the frame- 
work of hydrodynamic theory. However, our predictions 
are in disagreement in the BEC regime with the findings 
of Stringari 0, Heiselberg [|lj], and Manini et al [23, 
which predict a maximum in the transverse breathing 
collective mode frequency on the BEC side of resonance. 
This difference is probably due to our neglect of beyond 
mean field corrections |38l |39j . 

In conclusion, we have used a sum rule approach to 
study the collective oscillations of a two component 
trapped atomic Fermi gas in the BCS-BEC crossover 
region. This approach allowed us to simply compute 
the moments of the susceptibility without evaluating 
the full response function of the system. We explicitly 
treated the molecular bosons and were able to retain ef- 
fective interaction between Feshbach molecular bosons in 
the entire crossover region. Our calculation shows non- 
monotonic behavior of the breathing modes in the BCS- 
BEC regime. In the weak coupling BCS and BEC lim- 
its, our calculated breathing modes approach well known 
hydrodynamic results. At unitarity the breathing modes 
have same frequencies as those of a non interacting Fermi 
gas, as required by universality. The quadrupole mode 
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frequencies are independent of the interactions and the 
aspect ratio of the trap. Our predictions are in reason- 
able agreement with experiments and other theoretical 
approaches. 
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